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Abstract 

1 ^ I Following earlier studies, several new features of singular perturba- 

' tion theory for one-dimensional quantum anharmonic oscillators are 

computed by exact WKB analysis; former results are thus validated. 



X 



February 5, 2008 



This note continues our study [1] of singular perturbation theory in one- 
dimensional (ID) quantum mechanics using exact WKB analysis. Our focus 
', remains the v ^ 1 regime for the potentials V{q) = + vq^^ on the real 

line, with N > M positive even integers. Among those, the quartic os- 
cillator q'^ + vq^ has been a prime model for the mathematics of quantum 
perturbation theory [2, 3, 4, 5, 6]. Kawai and Takei [7] pioneered the use 
of exact WKB analysis in the latter context, followed by [8] (see also [9, In- 
troduction to Part I, and Pham's contribution], [10], and references therein). 
In spite of those successes, present exact- WKB quantization conditions (for 
g'^ + f g^, say) fail to tend toward their harmonic-potential {vq^) counterparts 

-'^Also at: Institut de Mathematiques de Jussieu-Chevaleret, CNRS UMR 7586, Uni- 
versite Paris 7, F-75251 Paris CEDEX 05, France. 
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as V ^ +00, be it analytically or numerically [11]. This worrying observa- 
tion triggered our present line of work (starting from [12, § 3]): to further 
probe how consistently exact WKB theory handles the perturbative {v ^ 1) 
regime. 

We are happy and honored to dedicate this work to Professor Kawai 
with gratitude, for his many essential contributions and leadership in exact 
WKB analysis, but also earher (with Professors Sato and Kashiwara) in hy- 
per/micro/function theory; this framework greatly inspired, and its Authors 
warmly encouraged, our first steps in exact WKB analysis [13]. 

Even though this work is thoroughly tied to [1] (with its bibliography), 
in § 1 we recall the main background and further strengthen the case for im- 
proper (divergent) action integrals like n(g) dq, where n(g) = {V{q) + A)^/^ 
is the classical momentum function. In § 2 we present new cases where 

n(g) dg can be computed exactly for some trinomial Il{q)^ '■ essentially 
the quartic case n(g)^ = + vq^ + A, for which Jq°° n(q') dg reduces to or- 
dinary (i.e., convergent) complete elliptic integrals. In § 3 we extend the 
main outcome of [1], namely the v +00 asymptotic expression of the spec- 
tral determinants D^(A,f) '= det^(— d^/dg^ + q^ + vq^^ + A) in terms of 

D^[K) =^ det^(— d^/dg^-|-g-^-|-A), to v — >• 00 in a complex sector. Thanks to 
this, finally in § 4 we demonstrate (currently provided 2M -|- 2 > TV) how the 
fundamental bilinear functional relation satisfied by does evolve into its 
counterpart for as — +00, in spite of a discontinuous jump at ^; = 00 
of the main parameter, the degree of the potential (from to M). 

1 Background (summarized) [1] 

Our model of quantum perturbation theory is the ID Schrodinger equation 



with N > M positive even integers. If we use the unitary equivalence (called 
Symanzik scaling) 



r d^ 



+ V{q) + a] *(g) = 0, g e R, V{q) = g"^ + vq 



(1) 



- dg2 



M 



,2/(M+2) 



r d2 



(Af+2)/(M+2) N 



(2) 



dg2 



+ uq + vq 



+ x"' uv 
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twice, atu = 1 and m = 0, the resulting right-hand sides imply that the opera- 
tor H = — d^/dg^+g^+fg^ is a singular perturbation of Hq = —d'^/dq^+vq^ 
for V ^ 1, and that the degree drops from N to M at v = +oo. 

From the classical dynamics we will use the momentum function (x i), 

n(g) = + A)^/^ (real in the classically forbidden region), (3) 

and its residue ReSg=oo n(g) = /9_i(0), a notation based on the expansion [12] 

(y(g) + A)-+V2^^^^(,)^P-iv. (^^oo); p = iTV, iTV - 1, . . . (4) 
p 

The spectrum of H is purely discrete, 0<Ao<Ai<---t+oo, and 
separates according to parity since V is an even function. Useful spectral 
functions (labeled by parity) are the generalized zeta functions, 



Z^s, A) Yl + (Re « > I + ^), (5) 

k 



odd 

1±/ 



and the spectral determinants D (A), defined through zeta regularization, 

logD^iX) = logdet±(^ + A) [-d,Z^{s,X)l^o, (6) 

where "s 0" implies analytical continuation in s. Scaling laws follow: 

det^[r(^ + A)] = r^*(*=°'^) det^(^ + A) (Vr > 0), 
where Z-(. = 0, A) . ± 1 

[11, equations (7), (30)] [12, equations (15), (27), (37)]. 

A more concrete realization of logD^ through (6) is, first to formally 
apply (d/dA)™ to (6) with the minimal m such that the result (oc Z(m, A)) 
converges, i.e., m > ^ + j^, then to integrate back: the separate knowledge 
that the A -\-oo expansion of logD^(A) shall only have "canonical" terms 
[14] [12, § 1.1.2] fixes the m integration constants. Here, > 4 implies 
m = 1 : specifically, 

AiogD±(A)^Z±(l,A) (8) 

converges according to (5), and logD^(A) is then the unique primitive of 
A) which is devoid of a constant (oc A°) term in its large-A expansion. 
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Classical analogs of those quantum determinants can be defined as well 
[11, 12, 10], through: logD^(A) = {the divergent part of logL'^(A) for 
A — > +oo}, or equivalently [12, § 1.2.1 and equation (46)] through: 

log(D+/D,l) (A) = logn(O) ^ ilogA, (9) 

/+00 f + OO 

Ii{q) dq = 2I, I / {V{q) + Xf' dg, (10) 
-oo -'0 

where this divergent "improper action integral" gets specified just like log D^: 
first, 

converges, then 7(A) is defined as that primitive of (11) which is devoid of a 
constant (oc A°) term in its large-A expansion. 

Improper actions as in (10) (i.e., along infinite paths) offer many benefits 
for asymptotic and exact WKB analysis. WKB solutions of (1) can now be 
defined intrinsically: e.g., as ^wkb(?) = n(g)~^/^exp J^^n(g')dg', unlike 
the traditional forms which awkwardly involve extraneous base points. The 
geometrical analysis no longer requires to set infinite paths apart as it used to 
[13, 17]. Moreover, the algebra itself is simplified; e.g., consider the full 

determinant D{\) =^ {D^ D^){\): previously, to get the large-A expansion of 
\ogD in the simplest case V{q) = q^^^ [13], we had to factor D = D^i a (a(A) 
is the "Jost function"), then expand logo using logo = J^^[U — Il]{q) dq 
where ^{q) — U{q)~^^'^ exp J U dq parametrizes an exact solution of (1), and 
finally obtain logDci by other means; now that the improper integrals (10) 
are allowed, all that condenses into a single identity (valid for general V): 

/ + 00 
U{q)dq. (12) 
-oo 

2 Explicit improper actions for trinomial n(g)^ 

In [1], we computed the improper action integral / = J^°° ^(q) dq in closed 
form for any binomial n(g)^ = uq^ + vq'^^; then (§ 4.2) we stated that we 
could no longer do so for a trinomial of the general (even) form Il{qy — 
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+ vq^ + A (with iV > M > 0), for which we just needed the v — > +00 
behavior of / anyway [1, equation (4.16)], reproduced as (36) below. 

It is nevertheless wrong to infer from the above that strictly no exact com- 
putations can be done in fully trinomial cases, and we now present several 
examples (still for even q^ -\-vq^ -\-\ with positive v, A). After recalling the 
closed-form results for binomials, we will quote another, trivial and degener- 
ate, instance: perfect-square trinomials. Then, our main new case will be the 
quartic anharmonic oscillator: we can reduce its improper action exactly to 
standard (i.e., convergent) action integrals, and therefrom to complete elliptic 
integrals [2, 15, 5, 16], as (30)-(32) below; we then verify the abovementioned 
large-i' behavior on this case {N — 4). Finally, the same approach must work 
for higher-degree polynomial V{q), converting /q^°° n(g) dq exactly into con- 
vergent hyperelliptic integrals (as studied in [17, 18]); but since the latter 
remain not so explicitly understood, we will skip this case {N > 4) here. 



2.1 Binomial Il{q)^ : exact evaluation 

For n(g)^ = uq'^+vq^^ (with > M > 0), n(Q') was exactly computed 
in [1, § 4.1]. We recall the main formulae for later convenience: 



(13) 
(14) 

(15) 



{uq'' + vq''Y/''dq''^]imIo{s), 

s~>0 

r+oo 

Io(s)= {uq^' + vq^'y/'-^dq (Re s > i + ^) 

Jo 

(iV-M) r(s- 1/2) 
where "s ^ 0" implies analytical continuation in s, with the result: 

'^W^^y^)''''"''"'^ (i=2(Ml)>l) when finite; (16) 
otherwise, i.e., when j = 1, 2, . . ., a further "canonical" renormalization yields 



I = 



2j/?-i(0) 



N + 2 



- i: - - ^(log2 + ,logK - )^ 



m=l 



/3-i(0) = (-1) 



i-i. 



(2j 



m=l 



(17) 
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We repeat from [1] the examples we will mostly need: 

r°°(wg^ + Xfl^dq ^ + ^) - ^) (^_^2) (18) 

Jo 2V7r 

- -| w-^/2A(log A - 1) (^ = 2) (19) 

/ (g^ + V)'/'dg = -i^;3/2, (20) 
io 

all based on (16) except (19), which uses (17) with j = 1. 

2.2 Perfect-square trinomials: n(g)^ = (g^^+\/A)^ (M even) 
This degenerate case trivially reduces to a binomial formula like (16), using 

for M > L > and w > 0. The singular formula (17) is never needed in our 
setting (n(g)^ even): for M and L even, no pole can appear in the numerator 
of (21) at s = ; but one appears in the denominator instead, leading to 

poo 

1= [(g^/2 + VA)2]V' dq = (for even N/2 > 0). (22) 
Jo 

2.3 The general even quartic case: Il{qy = + vq^ + A 

At present, we mean to exploit the large toolbox of results readily available 
for the complete elliptic integrals [19, 20, 21]: specifically here. 



K{k) = - t^){l - kH^)]-'^^dt, E{k) f 

Jo Jo 



^ ^ ^ ' dt, (23) 



(1 - t^) 

as functions of the modulus k ; the complementary modulus is k' =^ Vl — A;^. 
Main needed formulae 

- special values: [19, formulae 13. 8(5), (6), (15), (16)] 

K{Q) = E{Q) = ivr (24) 
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- derivatives: [20, formulae 710.00, 710.02][21, formulae 8.123(2), (4)] 

dK E(k) Kik) /dE E(k)-K(k). \ 

- expansions for k ^ 1' k' 0+ (implying E{1) = 1): [22, p. 93-94] 
[2, footnote 11 p. 184] [20, formulae 900.05, 900.07] 



K(k) = log ^ + ^ [log ^ - l] k'' + 0{k" log k') 



1 


[i 4 


1 




U 4 


13- 
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~ 2. 




~ 12. 



(27) 



selected transformation formulae: [19, Table 4 p. 319] 



1 + k' 



for k — 



k' 



1 (28) 



K(k)^k'K(k), E(k)^ 



m 

k' 



1 + k 

for A; = ^ . (29) 
k' 



Our closed- form result 

For non-negative v and A (as in [1], and mainly for simplicity), we find: 



{q' + vq^ + Xf'^dq (30) 

Jo 

{v>2y/X) : ^ \{v + 2V\Y/\2VXK{k) - vE{k)i A; = (^^^|^) (31) 

{v < 2x/A) : = i AV4[(2VA + v)K{k) - 2vE{k)l k = ^'^^'y^"^ • (32) 

Derivation. We first specify d//dA by means of (11) for V{q) — q^ -\- vq^. In 
contrast to (30), here the integrand {V{q) + A)~^/^ is integrable ed, q — oo 
in C, allowing to deform the path (— oo, +00) to a bounded contour in the 
complex g-plane: 



(33) 



where C is, e.g., a positive contour encirchng the pair of roots \q± of ^{qf 
{turning points) that he in the upper half-plane. We now prefer to pursue 
explicitly with v > 2\/A (and analytically continue the result to f < 2a/A 
later): then < g_ < cf. Fig. 1(a). The last integral in (33), being taken 
over a bounded path, admits a closed-form primitive with respect to A, as 

7(A) = 1 f {q' + vq' + Xy/'dq 

c 

'V?^ + V - A)'/' dq = -| q+[vE{k) - 2q^K{k)] (34) 
where q± = [|(w ± Vv^ - and k = [1 - ql/qlY'^, k' = q_/q+ 

[19, formula 3.155(1) for u = h and (amplitude) A = 7r/2] [2, formula (4.22)^]. 

We cannot rush to conclude that 1 = 1: the former contour deformation 
is ill-justified for the divergent integral / itself. On the other hand, we find 
that it simplifies future steps to use the transformation formula (28) which 
turns (34) into the expression (31), but still for I. 

Next, we continue (31) to the region {v < 2\/A} {k pure-imaginary) by 
means of the transformation (29), which results in the expression (32) again 
for I. Only then are we able to probe the A +oo behavior of / at fixed v : 
using k = ^ - 172^^'^^'^ + 0{X-^) and (25)-(26), we obtain 

/(A) ^ iA3/1X(^) + ^(^)-^.A-V2 + 0(A-)] 

+1 ^AV^[i^(^) + 0(A-V2)] - i ^AV^[^(^) + 0(A-V2)] 

-i^(7s)^'^'-i^^^'^' + ^(^~'^') (A-+00); (35) 

it has no constant (oc A°) term, hence indeed 1 = 1, the wanted canonical 
primitive as defined initially, cf. (10). QED. 

Remark 1. The trivial outcome (7 = I) seems to justify the above contour 
deformation directly for the divergent integral (30), but this is misleading: 
our n(g) kept a null residue /9_i(0), like all even T\.{q) with = mod 4 ; 
but generically, /9-i(0) 7^ (e.g., already for trinomial even n(g)2 but with 
N = 6, 10, . . .), and nontrivial integration constants /— / ^ ought to follow. 

^We think there should be no factor p^/-^ on the left-hand side of this formula. 
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Applications 

We can first verify (30)-(32) upon special cases, known earlier: 

- A = : ^ = -| v^/^E{l) = -i ^3/2 by (27) for E{1), cf. (20); 
-^; = 0: 7 = I ir(^) A=^/^ = by (25), cf. (18) for AT = 4; 

- V = 2VX : / = ^ v^/^[K{0) - E{0)] = by (24), cf. (22). 

But above all, we can use the exact expression (31) to check the v +oo 
behavior of / directly. Earlier, we predicted the asymptotic form for the 
general trinomial case to be, for v — > +oo at fixed A, [1, equation (4.16)] 

p+co r+oo r+oo 

/ (?^ + vq"" + A) dqr^ + vq"^) dg + / {vq^ + A) dg (36) 

Jo Jo 

+ 5m,2 ^^^Z^ v-'/'X{\og ^; + 2 log 2), 

where the first line is to be made explicit through (16)-(20), and 6 (last line) 
is the Kronecker delta symbol. 

However, our derivation of (36) was quite indirect, and lacked independent 
tests. Now the present results allow such a test: in the quartic case, we can 
directly expand I in its exact form (31) for v — > +oo, i.e., k — > 1~, and 

k' = 2{VX/vf\l + 2y/X/v)-^^^ ^ 0+. (37) 
Then, using (27), / = -| X^/'^k''^ [(2 - k'^) E{k) - k'^K{k)] expands as 

7 ~ -| A3/4 k'-i[l - I A;'2 - A(iog I - i) k'^ + 0{k'^ log k')] ; (38) 

the substitution of k' by (37) yields the desired v +oo expansion in terms 
of ti^/^^" and t>~-'^/2~" logf , n G N (no v^/'^Xogv term!). Remarkably, the 
next subleading term (of order v^/^) also cancels, so that finally 

7~ iA^;-V2(log(A/^;2) -41og2-l) [ + 0(^-^/2 log ^;)]. (39) 

This asymptotic equivalent then identically reproduces the prediction made 
by (36) for = 4 and M = 2 with the help of (19)-(20), which confirms our 
basic earher result [1, equation (4.16)]. 
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3 The V ^ oo behavior of the determinants 

We return to the spectral determinants of the quantum problem (1): 

D^{X, v) = det±(-dVd?2 + + vq^ + A), (40) 
which are entire functions of (A, w) e [24]. 

3.1 Review of the v — ^ +oo results 

Our key intermediate result in [1, equations (3.10-12)] was, for v — > +oo : 

D%iX, v) ~ e/o^°°(«"+-^"+^)^^^d« e- /o+-H-+A)V.d, det±(-d2 /dq' + vq^ + A). 

(41) 

Now the asymptotic formula (36) reduces this to 

D%{\ v) ~ e^('') qSmM\,v) ^g^±^_^2/^^2 ^ ^^M ^ ^^1^ (42) 

= /o+°°(g^ + ^g^)V2dg, ^(A,^) = ^(bg^ + 21og2) ^"V^a. 

(43) 

On the other hand, the exact scaling laws (2) and (7) for m = 0, plus 

/3-i(0) = 5m,2 a/2 for n(g) = {q^ + Kf'\ (44) 
entail (writing D%{K) = det^(-dVdg^ + q^^ + A)): 
det±(-dVdg2 + ^;g^ + A) = ^±i/[2(m+2)]^-5m.2.-V2a/8^± (^-2/(m+2)^). (45) 

our net asymptotic result was thus [1, equations (5.1-4)] 

^^(A, v) ~ e^(^) e^M,2Aix,v) ^;±i/[2(m+2)]^± ^ ^^^^ (4g) 

with I{v) = /+°^(g^ + vq^y/' dq given by (16) if j = ^ N, or by 

(17) otherwise, and 

Sm,2 A{X,v) ^ 5m,2 + log V + 4N log 2] A, (47) 

A = ^-V(iv^+2)A (= y-y^x in (47), used when M = 2). (48) 
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Figure 1: Plots of the Stokes geometry in the complex g-plane for n(g)^ = 
q'^ + vq^ + A and large complex v, ordered clockwise with increasing 6 — arg v 
{\v\ = 5, A = 0.5). The intermediate plots (b-e) set A = to emulate the 
|v| = oo regime at finite q; in that limit the Stokes curve S {bold line) stays 
linked to q = +oo (arrow) for 6* < (here = 27r/3, by (49)). 



3.2 Extension to a sector in the complex f-plane 

The key to our proof of the asymptotic formula (41) for positive v +oo was 
[1, § 3.2] that a solution "ifx^qjv) of (1) with a recessive WKB form for q 
+00 connects all the way down (in that WKB form) to a region {1 <^ q <^ 
^i/{N-M)y _ ^}^QTi-Q then tends to a similarly recessive solution \1/o,a(?) t^) of 
the uncoupled Schrodinger equation [— (d^/dg^) + vq^ + A]^o,a(q') = 0. 

In the complex domain, a simple sufficient condition for the WKB form 
to be preserved is for g e C to stay within one Stokes region of the momen- 
tum function Il{q) [13]. In terms of ^ = Bxgv, the above connection con- 
dition then becomes that the Stokes region containing {1 ^ e'^/(^+^)g <S 
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yi/{N-M)y (rotation given by the uncoupled equation) should link to q = +00. 
E.g., when V > the central Stokes region does include all of R, cf. Fig. 1(a). 
We now need to describe the Stokes geometry for Il{q)^ = + vq^ + A with 
complex V ^ 00; Fig. 1 illustrates the case N = 4, M = 2. 

When |w| — > 00 : the approximate factorization of n(g)^ as {q^~^ + v) x 
[q^ + \/v) makes M of its complex turning points qj shrink (x "in- 
ner" roots) and the other (A^ — M) grow (x v^/'y^-^) ^ "outer" roots); more- 
over, the central Stokes region contracts to a symmetrical pair of Stokes 
curves from g = for the zero-energy momentum 11^=0(5)^ = q^ + vq^^ , and 
we are to follow the (^-dependent) Stokes curve S which starts as = 1R+ 
when ^ = 0, cf. Fig. 1(b). In the large-f hmit, the connection condition 
is that 9 can be increased above as long as S remains linked to q — +00 
(Fig. l(c-d); the complex-conjugate picture results for ^ < 0). 

Following [17, § 3], the connection condition breaks (cf. Fig. 1(e)) when 
the action integral X — Jq° {q^ + vq^y^'^ dq becomes real, where go is the first 
outer turning point met hy S as 9 recedes from 0. That action, of instanton 
type [10], is computable in closed form: go = e~^'^/^^~^^ v^^^^~^\ and X = 
]^r(^)/r(^) [e-UM+2).^A.+2]i/[2(iv-M)]^ ^i^i^h turns real first 

at axgv — TT. Consequently, all of § 3.1 extends to the asymptotic sector 

E = {^;^oo, |arg^;| <e}, ® = (^9) 

Remark 2. Some examples (with A = 0) make us hope that our end asymp- 
totic formula (46) might actually hold up to | argv| < tt. 1) For Qi^(i') *= 
det^(— d^/dg^ + q^ + vq^), this was suggested by our numerical observations 
[12, equation (87) vs Fig. 1] that (^v^{v) behave analogously to the Airy 
functions Ai(f), Ai'(f) for v —00 as well (even though = 27r/3 only). 
2) The supersymmetric determinants det^(— d^/dg^ -|- g^ -|- vq^^'^~^) are 
known in closed form, essentially as inverse F-functions of v [12, equation (120)]: 
their large-^; asymptotics then amount to the Stirling formula, and the latter 
definitely holds for | arg^;| < tt (vs © = 7r/2). 



4 Asymptotics and the functional relation 

An early puzzle of general exact quantization conditions was their breakdown 
(both analytical and numerical) for potentials g^ + wg^ in the regime v 3> 1 
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(as seen for = 4 [11]). Naively, convergence to the elementary harmonic 
{vq^) behavior would have been expected. We can now show that singularity 
to be unessential: i.e., the original functional relations (Wronskian identi- 
ties) which produce those quantization conditions behave as well as possible 
when V — > +00 for any potential + vq'^ , currently under the restriction 
2M + 2 > N (which encompasses + vq'^, for instance). 



4.1 The basic Wronskian identity 

The spectral determinants for a general polynomial potential V{q) of de- 
gree N obey the bilinear functional relation: [11, equation (40)] 

where D"^^^^ are the determinants for the first conjugate problem: [24, § 7] 
V{q) ^ V^^\q) = e-'^^ V{e-'^^'^q) and A ^ = e"^^^ A, (51) 

with ifN = ' symmetry angle in degree N. (52) 

Equation (50) is but a Wronskian identity for the Schrodinger equation (1), 
yet it has a key dynamical role: while it seems underdetermined, it implies a 
complete set of exact quantization conditions, which then solve (1) exactly. 

A certain iterate of the transformation (51) is the identity, hence (50) has 
a cyclic symmetry group, specifically of order {^N + 1) when V is even. 

For the trinomial determinants (40), the first-conjugate parameters are 
^[1] ^ e-'^^ A, ^;W = e''^/^' v, with j = ^^j^^^^ as in (16). (53) 



4.2 The V ^ 00 transition 

According to (46)-(48) with A = -t;-2/(^+2)A , D%(X,v) for finite -y is a de- 
formation from D^j{A) at v = +00, but the key parameter in the dynamical 
functional relation (50), namely the degree of V, and often the residue /9_i(0) 
as well [1, § 3.1], sufi:er sharp jumps ai v = 00. It is then a non-trivial task 
to find out whether the basic identity (50) for continuously evolves into 
its counterpart for in the v — > -|-oo limit of (46), or not. 
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A AW 

Under (^^) i— (^^[i]) as in (53), the rescaled spectral parameter A maps to 

A ^ AW = exp[-^ i _ i (^^] A ^ e"''^- A ; (54) 

already this is the correct rotation angle for the limiting determinants D"^. 

To get the asymptotic form of (50) with = D'^(X, v), we let — > cxo in 
its left-hand side with argt; = — tt/ 2j, argw'^' = +7r/ 2j, and we invoke (46). 
The latter, by (49), requires n/ 2j < Q j > 1 or 2M + 2 > (otherwise 
the calculation will still work, but only formally until (49) extends to a wider 
sector). The left-hand side of (50) thus displays the asymptotic form 

exp [/(^;) + /(^;W)] exp (5m,2[^(A, i;) + A(AW, t;W)] x 

[z (AW)D- (A) - (A)D^(AW)], (55) 

We now evaluate all the terms in (55): first, 

I{v)<xv^ /(t;W) = -/(y)] if J ^ N (cf. (16), with /5_i(0) = 0) 

= /?_i(0) [logv + const.] if j e N (cf. (17), with /3_i(0) oc v^) 

^/(^) + /(t;W) = ^/5_,(0)if 

=> I{v) + 7(^;W) = i(^jv/3_i(0)/2 in all cases, cf. (53). (56) 
Next, just as for (54), 

^ ^ giViv/4 gi7r/[2(M+2)j] ^ ^i^N{N+2)/[4{M+2)] _ _ (-57-) 

Finally, and only relevant when M = 2, in which case A^^l = —A, 

^(A, v) + A(AW , ^W) = _ _N±2_ i ^ A = _ AL±2 i A = -i 7rA/4. (58) 

In the end, substituting (54)-(58) into (50) we indeed get 

e+i'^M/4^+(g-i<^MA)£)-(A) -e-i'^^/^D+(A)D^(e-''^^A) = 2ie^^'^'''^/\ 

(59) 

which is the correct form of (50) for D^j(A) = det^(— d^/dg^+g^+A) (whose 
/?_i(0) is given by (44)). For more details: if M > 2, see [23, equation (5.32)]; 
if M 2, then 

Dt{A) = 2^-^/2 V?/r(i±^), D^{A) = 2-V20F/r(3±A), (go) 
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and (59) with its "anomalous" right-hand side boils down to the reflection 
formula for the Gamma function; the harmonic-oscillator quantization con- 
dition can then also be recovered solely from (59) [12, Appendix A. 2. 3]. 

In conclusion, we have verified that the exact functional relation (50), 
governing both D^{X, v) and L>^(A) (cf. (59)), is compatible with the general 
perturbation formula (46), currently under the restriction 2M + 2 > N, or 
j > 1 (which includes the quartic oscillators): this further validates the exact- 
WKB description of perturbative regimes in [1]. Remaining desirable tasks 
are: 1) to lift the restriction j > 1 (e.g., by extending (46) to {| argf | < tt}, 
cf. Remark 2 in § 3.2); and 2) to find exact quantization conditions that 
themselves behave continuously in the zero-coupling limit (here, v — -|-oo). 
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